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The Kosterlitz-Thouless phase transition is described by 
the nonperturbative renormalization flow of the two dimen- 
sional </p*-model. The observation of essential scaling demon- 
strates that the flow equation incorporates nonperturbative 
effects which have previously found an alternative descrip- 
tion in terms of vortices. The duality between the linear and 
nonlinear cr-model gives a unified description of the long dis- 
tance behaviour for 0(A'')-models in arbitrary dimension d. 
We compute critical exponents in first order in the derivative 
expansion. 
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The Kosterlitz-Thouless phase transition may de- 
scribe the critical behaviour of various two dimensional 
systems. It poses a challenge to our theoretical under- 
standing due to several uncommon features. The low 
temperature phase exhibits a massless Goldstone boson 
like excitation despite the fact that the global U{1) sym- 
metry is not spontaneously broken by a standard or- 
der parameter. In this phase the critical exponents de- 
pend on the temperature. In the high temperature phase 
the approach to the transition is not governed by criti- 
cal exponents but rather by essential scaling. Modern 
nonperturbative renormalization group methods have de- 
scribed the characteristic properties of the low tempera- 
ture phase qualitatively Q . In this letter we extend this 
analysis to the essential scaling in the high temperature 
phase and provide for a quantitatively accurate discus- 
sion of the low temperature phase. We present a uni- 
versal renormalization group description of models with 
0{N) symmetry in arbritrary dimension. 

We employ the concept of the effective average action 
Tfc [§| 7 which equals the effective action F apart from the 
fact that in the former only fluctuations with momenta 
larger than k are included. For large momentum scales 
Ffc therefore equals some appropriately regularized mi- 
croscopic action while for fc = we have Fq = F. Thus F^ 
interpolates between microscopic and macroscopic scales. 
Since F is the generating functional for the IPI correla- 
tion functions it specifies directly the quantities of inter- 
est like the correlation length ^ = m . The flow of F^ 



is obeys an exact renormalization group equation 



5,F,[<^] = iTr((Ff [ 



Here F^ is the second functional derivative, dt denotes 
the logarithmic derivative k • d/dk and the trace (in mo- 



Rk) ^ ■ dtRk 



(1) 



mentum space) reads Tr — J2a=i I dl^q/ ■ The cut- 
off Rk{q^) suppresses the low momentum modes. We use 
a cutoff of the form 



Rk{q^) = Zkq^r{q^/k^) 



Zuq^ 



exj>{q'^ /k^ 



(2) 



where the wave function renormalization Zk will be fixed 
later. In order to solve equation numerically one has 
to truncate the most general form of F^. We introduce 
dimensionless, renormahzed fields, (pa = z\J'^k^'^~'^^^'^if!a^ 
p = {l/2)(pa<Pa and parametrize F^ in first order in a 
derivative expansion by 



Ffe = I d^'x^ k'^ukip) + ^k''-'zk{p)d^^adt,^a 



+ -^k''-'yk{p)d^pd,p + 0{d'') 



(3) 



The flow of Tk is then given by the flow of the quantities 
u, z and y which are dimensionless and depend only on 
the 0(iV)-invariant p and on k. We denote by k the 
running minimum of the potential Uk{p) and fix Zk by 
requiring Zk{n) = 1. 

The partial differential equations |^-^ 

dtu — —du + (d — 2 + vi)pu' 

+2vd{N ~ lX{w, z, ri) + 2vdli(,w, f , r;), (4) 

dtz — rjz + pz'{d — 2 + 7]) 

-{4vd/d)p^'^ {m^^Qiw, z,ij) 

~2mf -^{'w, w, z, z, rj) + rriQ ^iw, z, ry)} 

-2vdiz - z)p~'^{lf{w,S,r]) 

-{2/d){i-z)ll+\w,S,rj)} 

-2vdz'{{N ~l)lf{w,z,r]) 

-{8/d)nfj^{w,w,z,z,r]) 

+ {5 + 2z"p/z')lf{w,S,r]) 

-i4/d)z'plf+\w,w,zrz,v)}, (5) 

dfZ ~ i]z + pz'{d — 2 + 1]) 

-2vd{z' + 2pi")li{w,z,ri) 
+?>Vdpi' {iu" + 2pu"')4{w, z, f]) 
+ivd i2 + l/d)p{z')H^+^{w,2,r]) 
-{8/d)vdpiSu" + 2pu")'^rh'^{w, z, rj) 
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-(\Q/d)vdpi'{iu" + 2pw"')m4+^(w, rj) 

-{8/d)vdp{?Jfmi+\w,S,v) 

+{N-l)vd{ 

-{8/d)p{u"fmi{w,z,rj) 
- {16 /d)pu"z'm'l+'^{w, z, f]) 
-{8/d)p{z'fmi+\w,z,7j) 
+4(z - z)u"l^{w,z,r]) 

+4 {z'{~z -z) + {l/d)p{zr) lt^\w, z, ry)} (6) 

invoke the "threshold functions" ln^,n2^ ™ni.n2' "^ni,n2 
and nj^ „ defined by the integral 



X 



(7) 



with X ^ I, y{dypY, y{dypY, ydyp for I, m, m, n 
respectively. We have defined u' = du/dp and we use 
the shorthands v~'^ = 2'^+^'K'^/'^T{d/2), w = u' , w = 
u'+2pu", z = z+py, p{y) = y{z+r{y)), andp(y) = y{z+ 
r{y)). The derivative dt only acts on the fc-dependence of 
the cutoff i.e. dtp{y) = —y(rir{y) + 2ydyr{y)^ and we 
note that dtdyP — dydtp- Finally we abbreviate if^ q = if^ 
etc., where Iq is defined by the rule (p+w)^" — log(p+ 
w). The expression for the anomalous dimension r/ — 
—dt In Zk can be obtained from the identity dtZk{Kk) = 0. 
For = 1 one has y = Q, z — z. 

These equations are valid in arbitrary dimension. In 
the very simple approximation z — z — 1, u = iA(p — k)^ 
these equations already give a correct qualitative pic- 
ture for 0{N) symmetric models in arbitrary dimension 
[^. We show that the present version leads to quantita- 
tively accurate results. This is done by a numerical so- 
lution with initial values specified at a microscopic scale 
fc = A. For K ^ 1 the evolution is dominated by the 
iV — 1 Goldstone modes {N > 1). More precisely, the 
threshold functions at the minimum vanish rapidly for 
w = 2ku"{k) ^ 1. For N > 2 the coupling of the nonlin- 
ear cr- model for the Goldstone bosons is given by k~^. 

We concentrate first on c? = 2, where the universality 
of the /3-function for the nonlinear coupling implies for 
K ^ 1 an asymptotic form {N > 2) 



dtK = /3k 



iV-2 A^-2 



An 



167r2 



(8) 



In the linear description one can easily obtain an equa- 
tion for K by using dtUk{Kk) — together with equation 
(^. By evaluating the above equations for large k it is 
possible to compare with (|^) . Previously it was found ^] 
that in a much simpler truncation one already obtains 
the correct lowest order in the above expansion. In order 
to reproduce the exact two loop result one has, however, 
to go even beyond the truncation (||). 



TABLE I. Nonabelian nonlinear sigma model in d = 2. 
We show the ratio between the renormalized mass tur and 
the nonperturbative scale Aerge in comparison with the 
known ratio ^ invoking Ajj^: Cerge = ruR/ Aerge, 
^Ws ~ ^^/-^JTs' ~ Ttfl/fcs. We also display the ex- 
pansion coefficients for the beta function. 



N 


Cerge 




Cs 


/3i^' 




3 


2.81±0.30 


2.94 


1.00 


1.00 


0.79 


9 


1.22±0.03 


1.25 


1.05 


1.00 


0.84 


100 


1.08±0.04 


1.02 


1.06 


1.00 


0.87 



From equation (^) we expect that k will run only 
marginally at large n. As a consequence the flow of 
the action follows a single trajectory for large —t and 
can be characterized by a single scale. Notice that the 
perturbative /3-function vanishes for N = 2 since the 
Goldstone bosons are not interacting in the abelian case. 
Thus for large k one expects a line of fixed points which 
can be parameterized by k. This fact plays a major 
role in the discussion of the Kosterlitz-Thouless tran- 
sition below. It is responsible for the temperature de- 
pendence of the critical exponents. We have evaluated 
/3^=(3i^\N-2)/{ATT)+pi^^K-\N-2)/{16TT^) + ... nu- 
merically from the solution of eq. (|4|) - (||) and extracted 
the expansion coefficients for large n (see table |). 

For the nonabelian nonlinear sigma model in d = 2, 
N > 2 there exists an exact expression [|| for the ratio 
of the renormalized mass mn and the scale A-p-j which 
characterizes the two loop running coupling in the MS 
scheme by dimensional transmutation. The flow equation 
(Q) together with a choice of the cutoff Rk and the initial 
conditions also defines a renormalization scheme. The 
corresponding parameter Aerge specifies the two loop 
perturbative value of the running coupling similar 
to A-jg^ in the MS scheme. The numerical solution of 
the flow equation permits us to compute mn/ Aerge- 
(Two loop accuracy would be needed for a quantitative 




FIG. 1. The beta function for d = 2, N ^ 2. 
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FIG. 2. Essential scaling for d = 2, = 2. The renormal- 
ized mass uir is plotted as a function of ka = ka*—H{T—Tc). 

determination of Kerge/ ■^'jjg-) In table | we compare 
our results with the exact value of mn/ Kjjg. We also 
report the ratio mji/ks with kg defined by K{ks) = 0. 

The abelian case, = 2, is known to exhibit a spe- 
cial kind of phase transition which is usually described 
in terms of vortices . The characteristics of this tran- 
sition are a massive high temperature phase and a low 
temperature phase with divergent correlation length but 
zero magnetization. The anomalous dimension r] depends 
on T below Tc and is zero above. It takes the exact value 
77* = 0.25 at the transition. The most distinguishing fea- 
ture is essential scaling for the temperature dependence 
of rrifi just above T^, 



c = 



1 



(9) 



Our aim is to draw a picture of this transition from the 
flow equations without ever examinig vortex configura- 
tions. 

We have already mentioned the existence of a line of 
fixed points for large values of k, which is relevant for the 
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FIG. 3. Details of the beta function for d = 2, iV = 2. The 
curve is a fit to eq. (p]). 
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FIG. 4. Temperature dependence of the anomalous dimen- 
sion 77 for the low temperature phase, d = 2, N — 2. The line 
of fixed points is characterized by k and ends in the critical 
point for the Kosterlitz-Thouless phase transition. We also 
show the fiow towards the line of fixed points and the flow 
in the high temperature phase away from the critical point 
(left). The spacing between the points indicates the speed of 
the flow. 



low temperature phase. The contribution of a massless 
(Goldstone) boson in the renormalization group equation 
{w{k) — 0) is responsible for the finite value of rj. This 
in turn drives the expectation value of the unrenormal- 
ized field to zero (even for a nonvanishing renormalized 
expectation value k), 



po = njZk ^ K[klt.Y>. 



(10) 



We can observe this line of fixed points to a good approx- 
imation (cf. fig. 0), although the vanishing of /9k is not 
exact (we find -4 • Vd'^KT^ + The line of 

fixed points ends at a phase transition which corresponds 
to a microscopic parameter ka* ^ an initial value of the 
flow. In order to verify essential scaling we have to ex- 
amine the flow for values of ka just below that point, 
= KA* + <5'«A, <5ka ~ —{T — Tc). Then crosses 
zero at the scale ks and we find the mass by continuing 
the flow in the symmetric regime (minimum at k ~ 0). 
In figure ^ we plot (ln(TOfl/A))^^ against ka and find 
excellent agreement with the straight line (H). 

How does (3k, have to look like in order to yield essential 
scaling? Since there is only one independent scale near 
the transition, one expects mjj(T) = CgkaiT)^ where kg 
denotes the scale at which k vanishes, i.e. K{ks,T) = 0. 
For K close to and below we parameterize (this 
approximation is not valid for very small k) 



(11) 



For conventional scaling one expects C = and the corre- 
lation length exponent is given by v. Integrating equation 
( prj ) yields for C 7^ 0, (5k = k — k,: 
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TABLE II. Critical exponents v and 77 for d = 2. We 
compare each value with the exact result. 



iV 



0.70 
0.92 



0.75 
1 



0.222 
0.295 
0.287 



0.2083., 

0.25 

0.25 



TABLE III. Critical exponents v and r\ for d = 3 (see J^] 
for N = 1). For comparison we list in the third and fifth 
column an "average value" from various other methods |8|. 



N 


V 











0.590 


0.5878 


0.039 


0.292 


1 


0.6307 


0.6308 


0.0467 


0.0356 


2 


0.666 


0.6714 


0.049 


0.0385 


3 


0.704 


0.7102 


0.049 


0.0380 


4 


0.739 


0.7474 


0.047 


0.0363 


10 


0.881 


0.886 


0.028 


0.025 


100 


0.990 


0.989 


0.0030 


0.003 



ln(fc/A) 



(12) 



For k = kg the first term ~ (— is small and indepen- 
dent of T (since —5K,{ks) = k*) and equation ( p^ ) yields 
the essential scaling relation (^) for C — 1/2. Usually, 
the microscopic theory is such that one does not start 
immediately in the vicinity of the critical point and the 
approximation (|ll|) is not valid for fc « A. However, if 
one is near the critical temperature the trajectories will 
stay close to the critical one, Kc{t), with Kc(0) = ka*. 
This critical trajectory converges rapidly to its asymp- 
totic value At, and gets close to eq. ( pi] ) at some scale 
A' < A. As a result, one may use equation (|l^) only 
in its range of validity (fc < A') and observe that (5ka' 
is also proportional to Tc — T. The numerical verifica- 
tion of ( |lT|) is quite satisfactory: Fitting our data yields 
= 0.248, C = 0.502 and ly'^ = 2.54. The uncertainty 
for C is approximately ±0.05. The numerical values of 
and the approximation ( ITTI) are shown in figured. 

One can use the information from figure or ^n order 
to determine and therefore r/, = //(k,), the anomalous 
dimension at the transition. We plot rj against k in figure 
^. One reads off 77, = 0.287 ± 0.007 where the error 
reflects the two methods used to compute and does not 
include the truncation error. For k\ > ka* ot T < Tc 
the running of k essentially stops after a short "initial 
running" towards the line of fixed points. One can infer 
from figure ^ the temperature dependence of the critical 
exponent 77 for the low temperature phase. In summary 
all the relevant characteristic features of the Kosterlitz- 
Thouless transition are visible within our approach. 

We end this letter by reporting the values of the criti- 
cal exponents obtained in our approximation (^ - (^ for 
the "standard" second order phase transitions for d = 2, 
iV = 0, 1 and d = 3, TV > 0. In tables |l[ |!| they 



TABLE IV. 


Couplings 


for the scaling solution for d — 2 


and iV = 0, 1. 








N K, 


A. 




('^.) 


0.151 


5.33 61.6 


0.085 


1 0.265 


5.8 


B 65.4 0. 


868 


TABLE V. Couplings for the scaling solution for d = 3 and 


various A^. 








N 


A, 


7*3, 2l(«*) 


5.(k,) 


0.03009 


7.399 


78.84 0.192 




2 0.05984 


6.769 


51.25 -0.0415 


1.0602 


3 0.07651 


6.256 


39.46 -0.0920 


1.0695 


4 0.09414 


5.752 


30.52 -0.1107 


1.0789 


10 0.2162 


3.365 


8.17 -0.0584 


1.1144 


100 2.2313 


0.3779 


0.0947 -0.000759 


1.1468 



are compared with exact results or "averages" (only for 
the simplicity of the display!) of results from various 
other methods The agreement is very satisfactory! 
We also characterize in tables |^ the scaling solu- 
tion relevant for the second order transition by quoting 
K,, A, = ■u','(k,), M3, = u',"(k*) as well as z^k,) and 
z, (k,). We conclude that the first order in the derivative 
expansion of the exact flow equation for the effective av- 
erage action gives a quantitatively accurate picture of all 
phase transitions of scalar models in the 0{N) universal- 
ity class for arbitrary dimension 2 < d < 4. 
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